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1. Introduction 

The capability to grow nanoscale layers of atoms or molecules on a substrate forms 
an important part of modern nanoteclinology. Particles emitted from a source are 
deposited on the growing surface whose characteristics are determined by parameters 
such as growth temperature, surface diffusion, bulk relaxation, and particle attachment 
vs detachment rates. In addition to being clearly relevant for device applications, such 
surface growth problems also constitute an important class of generic nonequilibrium 
phenomena |2] which cannot be described in terms of standard equilibrium ensemble 
theory. 

A key goal of both theoretical and experimental investigations is the 
characterization and understanding of the resulting surface morphology, especially in 
terms of its statistical properties such as spatial and temporal height-height correlations. 
While local features may play an important role for device performance, large-scale, 
long-time characteristics offer better opportunities for theoretical understanding. They 
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are often universal, i.e., independent of microscopic details, depending only on global 
constraints such as symmetries and conservation laws. Here, we will explore the 
consequences of a reduced symmetry, which results if the particles arrive under an oblique 
(non- normal) angle. 

Let us briefly review different classes of models for surface growth. The first 
issue is whether surface overhangs and shadowing effects can be neglected. For models 
with near-grazing incidence, this is usually a poor approximation. Focusing on vapor- 
deposited thin films, Meakin and Krug IH Ej considered ballistic deposition under 
these conditions. They found that columnar patterns form which shield parts of the 
growing surface from incoming particles and characterized these structures in terms of 
anisotropic scaling exponents, differentiating parallel and transverse directions jSJ IH]. 
In contrast, for normal or near- normal incidence, overhangs and shadowing can often 
be neglected, and a description in term of a stochastic equation of motion for a single- 
valued height variable becomes possible. A second important issue is whether particles, 
once attached to the surface, can desorb again. If desorption does occur, the surface 
relaxation is not conserved, and the leading nonlinearity in the Langevin equation is of 
the Kardar-Parisi-Zhang form [Zj. Models of this type are most appropriate for chemical 
vapor deposition (see, e.g., 131^1) and sputtering processes (TUj. Anisotropies have also 
been included by different authors |TTl |T21 ESI- In contrast, desorption can often be 
neglected in MBE-type processes which occur at low growth temperatures (see, e.g., 
|14p. In this case, the surface relaxation terms in the Langevin equation conserve the 
number of particles. For models of this type, the universal behavior depends on whether 
the incident beam is normal to the substrate or not. For normal incidence, the system 
remains isotropic on large scales while oblique incidence generates anisotropies which 
persist under coarse-graining. While the former case has been investigated extensively 
(see, e.g., [13 HH]), only partial results are available for the latter Noting that 

the reduced symmetry generates novel leading nonlinearities suggesting potentially 
distinct universal behavior, we consider surface growth under oblique incidence more 
systematically in the following. 

Our analysis begins with a model first suggested by Marsili et al. \17^ to describe 
MBE-type or ballistic deposition processes with oblique particle incidence. Based 
on a coarse-grained description in terms of a Langevin equation, the model adopts 
an idealized perspective [IHl dl 120], neglecting particle desorption and bulk defect 
formation. As a result, the surface relaxation becomes mass-conserving and can be 
written as the gradient of a current, in an appropriate reference frame. In contrast, 
the noise term in the Langevin equation is not mass-conserving, since it models the 
random arrivals of particles. Anisotropies emerge naturally in the Langevin equation, 
since the particle beam selects a preferred ("parallel") direction in the substrate plane. 
Extending the analysis of fZj, we explore the consequences of an effective anisotropic 
surface tension. Generated by the interplay of interatomic interactions and kinetic 
effects, such as Schwoebel barriers, it plays a central role for the long-time, long-distance 
properties of the model. Due to the anisotropy, it will generically take different values. 
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denoted ry and r^, in the parallel and transverse directions. Depending on which of these 
parameters vanishes first, ripple-like surface instabilities are expected, aligned transverse 
to the soft direction. This leads to four different regimes with potentially scale-invariant 
behaviors. We analyze these four regimes, identify the scale-invariant ones, and discuss 
the associated anisotropic roughness exponents, using techniques from renormalized field 
theory. The original theory of J7] is recovered only if both couplings, r|| and Tj_, vanish 
simultaneously. 

The paper is organized as follows. We first present the underlying Langevin 
equation for a single-valued height field and briefly review the physical origin of its 
constituents. We also give a full discussion of its symmetries which are important 
for the subsequent analysis. We then turn to the renormalization group (RG) analysis. 
Identifying a set of effective coupling constants, and exploiting a Ward identity following 
from a (continuous) tilt invariance, we compute the naively divergent one-particle 
irreducible (IPI) vertex functions to one- loop order near the upper critical dimension. 
We analyze the four different cases separately: (o) both ry and r± are positive; (i) 
Ty remains positive while t± vanishes; (ii) ry vanishes while r_|_ remains positive; and 
finally, (iii) both ry and t± vanish simultaneously. A renormalization group equation 
allows us to derive the scaling properties of correlation and response functions for those 
cases which possess infrared (IR) stable fixed points. We identify a complete set of 
critical exponents, including four different roughness exponents for each case. Two 
of these characterize real-space scans along and transverse to the beam direction, and 
the remaining two are needed to describe scattering (i.e. momentum space) data with 
parallel or transverse momentum transfer. 

2. The model 

If surface overhangs and shadowing effects are neglected, the surface can be described by 
a single- valued height field, ^(r, t). Here, r denotes a d-dimensional vector in a reference 
(substrate) plane, and t denotes time, h is measured along the z-axis which is normal to 
the substrate plane. The no-overhang assumption can be justified a posteriori provided 
the interface roughness exponents are less than unity J7]. The time evolution of the 
interface is described by a Langevin equation, combining deterministic terms, G[h], and 
the effects of randomness, t], in the form 

dth = G[h] + r] (1) 

The conserved nature of the (deterministic) surface relaxation can be captured by 
writing G[h] = \/-F[h]. Several different terms, arising from surface diffusion and the 
incident flux, contribute to G[h] and are discussed in the following. 

The incident particle current has a normal component Jz and a component parallel 
to the substrate plane, Jy. For particles of finite size Tq, the effective flux must be 
measured at a distance Tq normal to the surface |2II. Since we focus on the long time, 
large distance characteristics of the growing surface, the corresponding contributions to 
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G[h] can be written as a gradient expansion, neglecting higher order terms: 

GfinAh] = -Jz + J|| ■ V/i + rj.y^h - r„J|| ■ {Vh) V^h (2) 

The first two terms of Eq. Q can be removed by transforming into an appropriate co- 
moving frame, via h{r, t) h{r + 3\\t, t) — Jzt. Next, we consider the effects of diffusion 
along the surface. This restriction generates a quartic term, of the form |22] 

Gdiff =-/iV^/i (3) 

Several comments are in order. First, we note the presence of the term ToJ^V'^h = rV'^h 
in Eq. Q. Similar terms can arise from a step edge (Schwoebel) barrier |2S1I211CS1, 
coming with a negative sign, or a surface tension with a positive sign. This allows for 
the possibility that the net contribution rV^/i might vanish, playing the role of a critical 
parameter. Even if we follow jT2j and set it to zero in the bare theory, it is actually 
generated under RG transformations and therefore intrinsically present. Second, the 
full rotational symmetry within the (i-dimensional space of the substrate is broken by 
the nonlinear term in Eq. As a consequence, any coarse- graining of the microscopic 
theory should give rise to different couplings in the parallel and transverse subspaces, 
and this is in fact borne out under the renormalization group. If the full anisotropy is 
incorporated into the theory, restricting rotational invariance to the {d — 1) -dimensional 
transverse subspace, the Langevin equation takes the form 

+ (Ai,9j/i + AxVi/i)+r7 (4) 

after some minor renamings. Here, d\\ (V_l) denotes derivatives (gradients) in the parallel 
(transverse) subspaces. It should be noted that the nonlinearity has also been split into 
two distinct terms, with couplings Ay and A_l. The Langevin noise ri{r,t) models the 
random particle aggregation at the surface, with zero average and second moment 

(r/(r, tMr', t')) = 2r'Sir - r') 5(t - t') . (5) 

In principle, one could also incorporate a noise term of the form —2{D\\d'^^ + D j_'V']_)6{r — 
r') 6{t — t') generated by (conserved) particle diffusion. In comparison to the aggregation 
noise, this contribution turns out to be irrelevant. However, one could imagine situations 
where diffusion is extremely fast, and aggregation occurs only very infrequently, so that 
Eq. (jH)) should be replaced by its diffusive counterpart. This "conserved" ideal MBE 
[23 1201 is not considered here, and will be left for future study. 

Eq. can be written in terms of surface currents j\\ and j_L, according to 

'y-^dth = - V±jx + ri (6) 

with 

jii = -9,1 (r||/. - /iii^j/i - /ix v» - ^{d\\h)' + ^{y^hf , 

j± = -V^(r^/i - /i^Vi/i - /xx9j/i) - A^(V^/i) {d\\h) . 

Eq. (jH), or its equivalent form, Eq. ©, is the starting point for the field-theoretic 
analysis. Its behavior is determined by the leading terms in the gradient expansion. We 
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Figure 1. Sketch of the phase diagram in {t^^,t±) space. EW denotes the 
Edwards- Wilkinson region. The hne labeUed 1 marks a hne of (possibly) discontinuous 
transitions, and line 2 corresponds to continuous transitions, both to inhomogeneous 
structures. The dashed lines denote very qualitative transition linesn between different 
structures, and mcr labels the multicritical point. 



can gain some qualitative insight by considering the hnear terms first. To ensure stabihty 
of the homogeneous phase, the coupling constants /xy and fi± must be positive, and the 
cross-coupling /Xx may not exceed —{lJ'\\fi±Y^^ ■ The two couplings r|[ and rj_ play the 
role of critical parameters. If both are positive, the nonlinearities become irrelevant, and 
the problem reduces to the (Gaussian) Edwards- Wilkinson equation [21] • In contrast, 
if one or both of them vanish, the surface develops characteristic spatial structures. 
These take the form of ripples (reminiscent of corrugated roofing), if only one of the two 
couplings goes through zero; if both couplings become negative, the surface develops 
mounds or "wedding cakes" p. Focusing only on the onset of these instabilities, four 
different cases emerge which are discussed systematically in the following section: (o) 
the linear theory ry > 0, r_L > 0; (i) a line of continuous transitions ry > 0, rj_ ^ 0; (ii) a 
line of possibly first order transitions t\\ —>■ 0, > 0; and (iii) the multicritical (critical 
end-) point t\\ — 0, r_L — > 0. The qualitative phase diagram is shown in Fig. (jT)). 



3. Renormalization group analysis 



3.1. The dynamic functional 

To unify the discussion, we first recast the Langevin equation (jH)) as a dynamic field 
theory. In this formulation, all statistical averages are expressed as path integrals, with 
weight exp(— 17). In order to treat correlation and response functions in an analogous 
manner, it is convenient to introduce a response field h{r, t), and formulate the dynamic 
response functional, following standard methods |2Zl I2H1 123 : 



J[h, /i] = 7 / d^'xdt^ h -i'^dth + d\\j\\ + VJ 



(7) 
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Before we turn to any computational details, we argue that the model is 
renormalizable. In other words, in each loop order we need to add only those primitive 
counterterms which can be absorbed in a renormalization of the fields and the coupling 
constants. As a first step in this procedure, we consider the symmetries exhibited by 
((Zj) since they place constraints on possible counterterms. The invariance against a shift 
of coordinates in the z-direction leads to the symmetry h{r,t) — > h{r,t) + a. Hence, 
the surface height h arises in the functional only in form of derivatives. Second, the 
functional is invariant under tilts of the surface by an infinitesimal "angle" b, i.e., 
h{r, t) — » /i(r, t) + b ■ r provided the tilt is accompanied by an appropriate change of the 
parameters, namely, t\\ r|| = t\\ — h\\\\\ and t_l ^ = r_L — 6||A_l. Further, particle 
conservation on the surface lead to the invariance transformation h{r,t) h{r,t) + c, 
h{r, t) —>■ h{r, t) + 2c7t . Therefore, we have a (d + 2)-dimensional (infinitesimal) 
symmetry group 

/i(r, t) h{r, t) + a + h-r + 2c'yt , h{r, t) h{r, t) + c , 

which does not change the form of the dynamic functional or the functional integration 
measure, provided it is accompanied by an appropriate transformation of the parameters 
r|| and t±. In addition, we have a discrete symmetry, namely, inversion invariance in 
the {x\\,z) subspace: 

h{x\\, r_|_, t) —h{—x\\, r_|_, t) , h{x\\, r±, t) — > — xy, r_L, t) . 

Finally, of course, we have isotropy in the transverse subspace. These symmetries can 
be considered as the fundamental defining elements of ideal MBE processes. 

The renormalizability of our model can now be argued, following standard methods 
[SHI EII- In each loop order, the successive construction of the perturbation series 
produces only primitive ultraviolet (UV) divergences, provided that all counterterms 
constructed in lower orders are included to cancel non-primitive divergences which 
appear in subdiagrams. If the regularization procedure respects the symmetries of the 
model, then the remaining primitive divergences also preserve the symmetries. Thus, 
these divergences have the same form as the relevant terms of the initial model and 
their counterterms can be absorbed in a renormalization of the initial parameters. 
This presumes, of course, that all relevant interaction terms (i.e., composite fields with 
coupling constants of positive naive dimension) with the correct symmetries have been 
included in the initial dynamic functional, and all irrelevant couplings have been omitted. 
With regards to our model, a straightforward analysis of the relevant diagrams shows 
that we have indeed captured all relevant terms in the original Langevin equation Q. 
Hence, the model is renormalizable. 

3.2. Elements of perturbation theory: the four different cases 

In this section, we assemble the basic components of the perturbative analysis, leaving 
technical details to the Appendix. We first write the dynamic functional, J^[h, h] as the 
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sum of two parts. J^o[h, h] contains the Gaussian terms of the theory, 
Jo[h,h]=^ J d''xdtJjii^-^dt + i^iidl + fi±i\/lf 



+ 2/Xx9jVi - riiaj - r^Vl)h - j (8) 

while J7i[/i, h] takes the nonhnear interactions into account: 

J^[h, h] = -j j (fxdt h {d\\h) {\\\dlh + x^ylh) (9) 

Arbitrary correlation and response functions can now be computed as functional 
averages with statistical weight J[h, h\. 

(o) = j oe-^I^''^!. 

Due to the anisotropy, it is possible to rescale parallel and transverse lengths 
independently. Considering a simple rescaling of parallel lengths first, via x\\ — > axy, 
the functional remains invariant provided h a^^^'^h, h a~^/'^h and /iy 
r|| cP"T\\^ A^x ct^/Wx while Ay a^/^Ay and A_l — > a^^'^\^_. Likewise, if only 
transverse lengths are rescaled, via r_L l3r±, the functional remains invariant provided 
h (3~^'^-^y^h, h and /ix ^ /3V±, ^± /5^^±, l^x ^ /S^x while 

A|| /^'•'^""^•'/^Ay and A_l — >■ P^'^^^^^'^X^. In addition, there is an overall length scale 

which accounts for the dimension of fields and coupling constants under the 
renormalization group. When we consider the three different cases below, these scale 
invariances will allow us to define appropriate effective expansion parameters, in order 
to eliminate redundant couplings |32j . 

Our goal is to compute the Green functions G'^^({r,t}) = ([/;,]^[/;,]^)™™, i.e., the 
cumulant averages of all possible products of fields h{r,t), h{r,t), or equivalently, the 
associated one-particle irreducible (IPI) vertex functions, r^^({r,t}), with N (N) 
h {h-) amputated legs, order by order in a diagrammatic perturbation expansion. 
The notation {r, t} is short-hand for the full space- and time-dependence of these 
functions. As a first step in this process, we collect the elements of perturbation 
theory and introduce their graphic representation. Using to denote averages taken 
with the (Gaussian) weight j7o, we will need the bare propagator, GqIQ;^) — q') = 
(/i(q, t)h{—q', 0))o , as well as the bare correlator, Co(q, t) (5(q— q') = {h{q, t)h{—q', 0))o , 
both in spatial Fourier space. In the following, we write ... as short-hand for 

I^It--- Defining 

r(q) = mi + /^±(qi)' + 2/ixgyqi + rygj + r^qi (10) 



we find 



Go{q,t)=9{t) exp(-r(q)7t) , 

Co(q,t) =r(q)-^Go(q,|t|). (11) 
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Here the Heaviside function 6{t) is defined with ^(0) = 0. Diagrammatically, these two 
functions are represented as hues: 



Go{q,t) = t. 



q 



Co(q,t) =t 3 

Turning to the interaction terms, it is convenient to rewrite them in a symmetrized 
form. In Fourier space, the expression for the three-point vertex reads 



^(qi, q2, qa) = h{h qiw q2\\ gay + 



+ A, 



qi± ■ (92|| q3± + q2± qz\\) - qi\\ q2± ■ q3± } (12) 



with qi + q2 + qs = 0. Our sign convention is such that all momenta attached to a 
vertex are incoming. Diagrammatically, the whole vertex is represented as 

^(qi,q2,q3) 

Returning to Eqs. (jH]) and ©, we now perform dimensional analysis. In the (infrared) 
limit of small momenta and frequencies, the Gaussian part of the dynamic functional is 
dominated by different terms, depending on the behavior of the control parameters ry 
and r^. 

Case (o). If both ry and rj_ are finite and positive, the theory turns out to be purely 
Gaussian. Quartic derivatives can be neglected in the infrared limit; as a result, r(q) 
simplifies to r*^^)(q) = ryg^ + T_|_q^. Since both parallel and transverse momenta 
appear only quadratically here, it is natural to choose a momentum scale k such 
that q\\ oc |q^| oc k, and ry oc oc Time 7t scales as k"^, and the fields have 
dimensions h{r,t) oc k('^~2)/2 ^^^^ h(r,t) oc k*^'^"''^^/^. The nonlinear couplings scale as 
A|| oc Ax oc and are therefore irrelevant in any dimension d > 0. The resulting 

theory is a simple anisotropic generalization of the Edwards- Wilkinson equation [2^1 , 

^-^dth = T\\df^h + Tj_Vlh + 7] . (13) 

The anisotropics in the quadratic terms affect only nonuniversal amplitudes and can be 
removed by a simple rescaling, without losing any information of interest. As is well 
known, the two-point correlation function scales as 

C(r,t) = |r|'-''c(t/|r|') (14) 

from which one immediately reads off the (isotropic) roughness exponent x = (2 — (i)/2 
and the dynamic exponent z = 2. For a detailed discussion see [2]. Since this case is so 
familiar, we will not consider it any further. 

Case (i). If ry remains finite and positive but tj_ is infinitesimal, the two leading terms 
in the dynamic functional are ^±_h{W\Yh and T\\hd?^h. Hence, even in the Gaussian 
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theory, parallel and transverse momenta already scale differently, reffected in the choice 
|q_L| oc K and q\\ oc h?. If we introduce an anisotropic scaling exponent A via q\\ oc k^"*"^, 
we recognize that A = 1 at the tree level. Continuing with the two leading terms, we 
note T\\ oc /i_L oc and set /ij, = 1 via a transverse rescaling with an appropriate (3. 
Time scales as and /iy oc /ix oc are both irrelevant. The only strongly relevant 
parameter of the theory is oc k^. Introducing the effective dimension D = ci + 1, 
one finds h(v,t) oc k^^~^^^'^ and h{v,t) oc k^^'^^^^'^. For the nonlinear couplings, we 
obtain A_l oc k^*^"^^/^ and Ay oc k^^/'^. Since D is clearly positive, the coupling Ay 
becomes irrelevant. The upper critical dimension dc for the theory is determined by A^, 
via = 4 — D which leads to = 3. The invariant dimensionless effective expansion 
parameter is t^^^^X±/k^'^~'^''^'^ as shown by the rescaling A_l a^/^A_L, ry — > a^ry. The 
resulting propagators are controlled by r(q) = T_Lq5_ + rygjj + (q5_)^ where the (q5_)^-term 
determines the UV behavior and t± plays the role of the IR cutoff. At the tree level, 
this case corresponds to a critical line parameterized by ry. 

Case (ii). Here, t± remains finite while ry vanishes. The Gaussian part of the functional 
is dominated by hd^^h and T±hV\h. Again, parallel and transverse momenta scale 
differently, but their roles are now reversed: gy oc k and |q_L| oc with A = 1 at 

the tree level. Both /xy and scale as while fi± and /ix are irrelevant. Defining 
the effective dimension as D = 2(d — 1) + 1, one still has h(r,t) oc k^^'^^^"^ and 
h{r,t) oc K*^-^"^^^/^; we also recover 'jt oc Now, however, the strongly relevant 

perturbation is ry oc k^. The two nonlinearities switch roles so that Ay oc and 
A_L oc K"^/^. Clearly, A_l is irrelevant in all dimensions while Ay becomes marginal at 
the upper critical dimension dc = 7/2. The invariant dimensionless effective expansion 
parameter follows from the rescalings as /Uy ''^'^r^*'"' ^^^^ X\\/ n^'^^'^'^^/^ . The momentum 
dependence of the propagators simplifies to r(q) = rygjj + T_i_q]_ + /iygy. As in the 
previous case, at the tree level this situation corresponds to a critical line parametrized 
by T±. However, we will see below that the order of the transition may well become first 
order once fluctuations are included. 

Case (iii). Finally, we consider the multicritical point where both t± and ry vanish. 
Both momenta scale identically, as gy oc |qj_| oc n, so that A = at the tree 
level. Again, we may set fi± = 1. One obtains '-ft oc ry oc t± oc k^, and 

h{r,t) oc K*^-^"'^)/^, h(r,t) oc k^^"*"^^/^, with D = d. The full propagators come into 
play. Both nonlinear couplings. Ay and A^, have the same upper critical dimension 
dc = 6. The effective expansion parameters are w = fix/^/J^, u\\ = /iy ''^^Ah/k'^^"'^)/^, 
and u± = /i|| ^^^X±/k^^~'^^^'^. In this case, the anisotropy exponent A vanishes at the tree 
level. 

In the following, we analyze the three nontrivial cases in a one-loop approximation, 
using dimensional regularization combined with minimal subtraction [201 IM]- The 
essential components are the IPI vertex functions r^^({q, u;}). Focusing on the 
ultraviolet singularities, only those F^ ^ with positive engineering dimension are to be 
considered. Taking into account the symmetries and the momentum-dependence carried 
by the derivatives on the external legs, the set of naively divergent vertex functions is 
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reduced to Fi i and Fi 2. Specifically, at the upper critical dimension of cases (ii) and 
(iii), Fi 1 is quadratically divergent and F12 is marginal. In case (i), Fi 1 is already 
marginal after accounting for the external momenta. 



3.3. Perturbation theory. 

In this section, we summarize the calculation of Fi 1 and Fi 2, to one loop order. Details 
are relegated to the Appendix. 

The main features which recur in all three cases are the following. If Fi 1 is 
quadratically divergent, this divergence is first removed by an additive renormalization, 
i.e., a shift of the critical parameter (s). Specifically, the true critical point is located 
through the singularity of the parallel and/or transverse static susceptibilities, x\\ 
X±, defined by 

7/^11 XiT^ = lim a|7^Fi,i(g||,q_L = 0,iu = 0) 

7Xl' = lim(qi)-iFi,i(g|| = 0,q^,u; = 0) (15) 

q±-*o 

We find that nontrivial shifts are only required for cases (ii) and (iii), and only for ry. 

The remaining logarithmic divergences are then computed using dimensional 
regularization, so that ultraviolet divergences appear as simple poles in e = dc — d. In a 
minimal subtraction scheme, we focus exclusively on these poles and their amplitudes. 
Since the nonlinearities are cubic in the field, the expansion is organized in powers of 
My, u]_, and u\\u±; i.e., the first correction to the tree level is always quadratic for Fi^i 
and cubic for Fi 2. The tilt invariance leads to a Ward identity connecting Fi 1 and Fi^2; 
namely. 



_d_ 
. d 



r f ' '\ \ 9Ti,i{q,uj) gFi,i(q,^) 
Fi,2(q,^;q,^) = Ay — + Ax (16) 

Fi,2(q,u;;qV) =0 

q',(^'=0 

Here, (q, a;) and (q u') denote the momenta and frequencies labelled with subscripts 
2 and 3 in Eq. p2j) . By virtue of this symmetry, it will be sufficient to compute Fi 1 in 
order to extract the renormalization of all parameters, including those of the coupling 
constants Ay and A_l. More specifically, the tilt transformation h{r,t) — > h{r,t) + b ■ r, 
''"II ^ ''"II ~ ^11-^11' ''"-L ^ ''"-L ~ ^11 -^-L shows that the parameter b renormalizes as the field h 
itself. Hence, the term Ay h renormalizes exactly like ry, and X±h like t±. 

Next, we collect a few general properties of the diagrams contributing to Fi i(q, cj) 
which follow from the special form of the interactions, Eq. 0. We first consider 
diagrams where the external h-leg carries two transverse derivatives. Such diagrams 
can arise only from the interaction Xj_h {d\\h)'V\h. Thanks to particle conservation, 
all these graphs sum in such a way as to generate at least one spatial derivative on 
the external h-leg. Moreover, either parallel inversion or transverse rotation symmetry 
guarantees that, in fact, there are at least two such derivatives. Thus, there are no 
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corrections to the term '~fT_i_h'V']_h in JTq, Eq. (jSI). Similarly, since external h- and h-legs 
always come with at least one spatial derivative, there are no perturbative contributions 
to the terms hdth and 7/1^. Hence, these three terms are not renormalized. Further, if 
the coupling constant is zero, as in case (ii), only diagrams with parallel derivatives 
on all legs are generated, and so operators carrying no parallel, or mixed parallel and 
transverse, derivatives suffer no renormalization. In contrast, when A_l 7^ 0, there is no 
such simplification since this vertex can contribute to purely "parallel" terms. If both 
coupling constants Ay and A_l are nonzero, they mix under the RG, and we have to 
renormalize the couplings by matrices, rather than simple multiplicative factors. 

With these considerations in mind, we propose the following renormalization 
scheme: 

h^h =Z^''^h h^h = Z~^l'^h 7^7 = ^7 

T^^f^ = Z'^T± T\\ f|| = Z''^{ZrT\\ + yUy^Vr^) + f||_c 

/ix ^ Ax = Z^'^Zy.fly, ^\\ = (17) 

Ty^c accounts for the shift of the critical point in cases (ii) and (iii). The Ward identities 
then lead to the renormalizations 

X±^X± = Z-'^l'^\^ A|| ^ A|| = Z-^l'^{Zr\\\ + Y\i_) . (18) 

The five renormalization factors Z, Z^-, Zy., Z^, and Y must be calculated from the 
e-poles of the dimensionally regularized perturbation series of ri i(q, cj). In minimal 
subtraction, they are of the form Zo = 1 + Yl'^=i -^o^V^") ^ = Xl^i ^'•"V^"; where the 
Z^"^ and are expansions in the renormalized effective coupling constants. Here, the 
subscript o stands as a placeholder for any member of the set of indices. The individual 
characteristics of the three cases of interest are now treated sequentially. 



Case (i): r_|_ ^ and ry 7^ 0. This is the simplest non-trivial case. Only one pa- 
rameter, r_L, needs to be tuned to access criticality. Since Ay is irrelevant, it may be set 
to zero. Neglecting all other irrelevant terms as well, the functional simplifies to 

J = 1 j c/'^xrft|/i[7-i9t+(Vi)2-r||9j-r^Vi]/i-/i2-A^/i(9||/i)Vi/i}(19) 

In this case /iy = /ix = Ay = F = 0. All divergences are parameterized hy e = — d 
where (ic = 3. 

Following standard methods (for details see [Appendix A ), we obtain the following 
expression for the singular part of the two-point vertex function, 

ri,i(q, u) = ^u + 7[r||g; + (qi)^ + r^qi] + 7^[2rygJ - (qi)^] + O{u'){20) 

where 

u = A,T-^^^X^K-'/^ (21) 

is the effective expansion parameter. A,, is a geometric factor which appears in all 
integrals: 

' {27Ty^ \ 2 ; 
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and Sd is the surface area of the d-dimensional unit sphere. 

The renormahzed vertex function Fi i is defined by demanding that 

ri,i(q,c^,7,^±,^lhM, k) = fi,i(q,u;,7,f^,f||, V) 

be pole-free. One finds 

Z = l + ^ + 0{u'), Zr = l-^ + 0{u'). (22) 

The corresponding Wilson functions are defined as the logarithmic derivatives of 
the associated Z-factors, at constant bare quantities, i.e., (o = t^d^^hi Zo\hare ■ The 
logarithmic derivatives of the control parameters and of the relaxation coefficient are 
then given by 

KS^lnrillbare = C " Cr , 

fi:9«;lnrj_|bare = C, (23) 

K9«,ln7|bare = "C, 

The flow of the dimensionless effective coupling constant u under renormalization is 
controlled by the Gell-Mann-Low function, 

/3 = fi:9«M|barc = M[-| + ^(C + Cr)]. (24) 

The renormalizability of the theory demands that all the Wilson functions be free of 
£-poles; moreover, in minimal subtraction it can be shown [HO] that the perturbative 
corrections do not contain e any more. It is easy to demonstrate that all Wilson functions 
are determined by the flrst term, zi^\ in the Laurent expansion of the Z-factors. So, 
we obtain 

m5Z« 4, 
an) =-2^ = -y + 0(-^)- 

/?(«.) =u[-'- + + 0(u*)] . (25) 

The renormalization group equation (RGE) for the Green functions simply states 
that the bare theory is independent of the external momentum scale k: 

= K^G^^^({r,t};fj_,f||,7,A_L) 

= «^Z(^~^)/2G^^^({r, r^, r,,; u, 7, n) . (26) 
Explicitly, this partial differential equation becomes 

It will lead to asymptotic scaling, provided I3{u) possesses an infrared stable fixed point 
M*, i.e., a solution of (3{u*) = with f3'{u*) > 0. To one loop order, Eq. ()25|1 allows for 
a single fixed point with the desired properties, namely, 

= ±4^1 (1 + 0(5)), (28) 
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where the sign is determined by the sign of the initial couphng constant Aj^. Clearly, 
this u* is non-zero in the e-expansion. Under the (reasonable) assumption that the 
full (resummed) series for u*{e) is non-zero at the integer values of e corresponding to 
physical dimensions, Eq. (j^^ gives us the exact relation 



au*)+Cr{un = f (29) 

In the following, Wilson ^-functions, evaluated at the fixed point, will simply be denoted 
by a superscript *. 

The solution of the RGE (j2Zj) is easily found by applying the method of 
characteristics. At the stable fixed point u* we obtain the scaling solution 

G^^^({r, t}; r||, 7, k) = /f*(^-^)/2G^^^({r, t}; Fr^, l<'-^'r^^j-C^^ 

where / is an arbitrary flow parameter. The parallel and transverse scale invariances, 
together with dimensional analysis, give us an additional scaling relation for the Green 
functions, namely 

G^,^({r,t};r^ ,r||,«^,7,K) = «(A^+^)/2^A^(^-4)/2+iV(D+4)/2 

Combining these two relations, choosing a = /K^"^*)/^ and recalling D = d + 1, we arrive 
at the final scaling form for the Green functions. Suppressing unneeded arguments, we 
can choose to write them in the form 

Gjv,^({r, t}; rx) = /'^-.-G^,^({/i+^X||, /r^, IH}; r^'^'T^) (30) 

where the scaling exponents A and rj can be expressed in terms of Wilson functions, 
evaluated at the fixed point: 

A = l + ^^^, r, = C (31) 

Since we have only two independent renormalizations, the remaining exponents follow 
from scaling laws, namely, 

z = 4-77, l/z/ = 2-r/, (32) 

and 

^N,N = y(rf + A- 4 + r^) + ^(rf + A + 4-r^) (33) 

Eq. ()29p. provided it holds, can be exploited to give another exponent identity: 

A + 7/ = 2 - rf/3 . (34) 

To conclude, only a single exponent, e.g., rj has to be computed order by order in 
perturbation theory. Then, all others follow from exponent identities which are exact, 
at least within perturbation theory. Our one-loop calculation results in 

T] = -2£/3 + 0{e^) . (35) 
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Case (ii): ry and t± ^ 0. This is the second non-trivial case. Neglecting irrelevant 
terms, the dynamic response functional simplifies to 

J[h,h]=-f j (fxdt^h 7-^9t+/i||9j-r||9j-rxVl h-h^-\\\h{d\\h)dlhY'iQ) 

In this case, we define e = dc — d with dc = 7/2. 

Since ri^i(q, tu) is quadratically divergent in this case, an additive renormalization 
of the critical parameter, t\\, is required (see [Appendix B ). Once this is accounted for. 



we obtain the following one- loop expression for the singular part of Fi i: 

ri,i(q, w)poic = + 7 (^/iiigf + rygj + r^ql) + 7^^ (^mfl + ^II^S) + • • • (37) 
where u is the effective expansion parameter now defined as 

„^5^^-7/8^;('^-i)/4A||«:-(^-2^)/^ (38) 
with a geometric factor 

j2 _ Sd-i I d — 1\ ^ I d 



Bt = ^r^V r - r 1 + e . 

Eq. ()37j) gives rise to renormalizations of /iy and t\\. Since the vertex Ay carries only 
parallel momenta in this case, there are no contributions of order q^^ at any order, 
implying the absence of field renormalizations. Hence, Z = 1. The remaining poles are 
absorbed into the renormalization factors 

^. = l-^ + 0(n^), 

Zr = l-- + 0{u') . (39) 

Defining the Gell-Mann-Low function of u, as well as the Wilson functions, in analogy 
to the previous case, we obtain 

C = l^df. In I bare = + 0{u'^) , 

= nd. In Z.lbare = + O(m^) , (40) 

7 1 

P = Ibare = u[ - E + -(^^ - Cr] = u[ - £ - -U^ + 0{u^)] 

o o 

The fixed point equation I3{u*) = does not yield a stable real fixed point, at least to 
this order in perturbation theory. It remains an open question whether such a fixed 
point, and the associated scaling properties, might emerge at higher orders. Restrict- 
ing ourselves to our current results, the absence of an infrared stable fixed point might 
suggest a first order transition. A more detailed analysis of the underlying mean-field 
theory or a careful computational study would be required to shed more light on this 
issue. Both are beyond the scope of this article. 

Case (iii): and t± 0. Finally, we turn to the analysis of the multicritical 

point. This situation was previously considered in ^Jj, using a momentum shell 
decimation scheme. This procedure requires considerable care for field theories with 
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strong anisotropy and nonlinearities carrying multiple derivatives, since the corrections 
depend on the way in which the hard momentum cutoff is implemented. Moreover, 
the hard momentum shell cutoff introduces long-ranged correlations on the scale of 
the cutoff momentum which must be handled very carefully. Whatever the source of 
the discrepancies, even after meticulous checks we were unable to reproduce the earlier 
results. 

The full functional, Eqs. (jHl) and Q, as well as the full (bare) propagators and 
correlators, Eq. (|TT|) . now come into play. To ensure the stability of the critical theory 
at the tree level, we demand fi\\ gj| + + (qi)^ > 0. This limits the physical 

range of /iy and fix to > and fix > ~y/JH- To complicate matters further, both 
nonlinear couplings. Ay and A^, are marginal at the upper critical dimension dc = 6. 
However we know that, thanks to the Ward identity, all renormalizations can still be 
obtained from the two-point function, Fi i. 

A first analysis of diagrams contributing to ri i(q, cj) shows that two of them are 
quadratically divergent. Since both of them carry a momentum prefactor of gj^, they can 
be absorbed in an additive renormalization, i.e.., a shift of the bare control parameter 
ry, as remarked before. No shift of the bare t± is needed. 

Leaving details to [Appendix C[ our one-loop result for Fi i is now considerably 
more complex, thanks to the presence of both vertices. In order to eliminate redundant 
parameters, it is again convenient to define invariant dimensionless coupling constants, 
guided by the parallel and transverse rescalings: 



nil 



w = 



1/2 



where the constant is defined by 

7rF(l + e/2) 
" {27iY 32V2 

The definition of the effective coupling constants contains a suitable w-dependent 
denominator, which is common to all integrals. Collecting the results from [Appendix CI 
we find the singular parts 



(41) 



Fi,i(q, w)poic = + 7 g||(/X||g|| + ry) + qx(qx + r^) + 2/ixg||q± 



1/2 2 

7/i|j q\\ 



Mil 



Bi -1/2 



Bi 



2 (B2 -1/2 



^B, 



2C1 , 2 '-'2 
Mm h M I 



e 



1/2 2 2 

7/^11 ?i|q± 



u 



,D,^ ^2 



Do 



I 2\2 2 ^1 

7(q±) u^— 
D3 



(42) 



U\\U_i_- 



e e e 

Here, Ai, Bi, D3 are w-dependent functions whose explicit forms are listed in 



[Appendix C[ 
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According to the general renormahzation scheme, Eq. (jl7|) . the full set of five 
renormalization factors now comes into play. At one loop order, we read off 



?/2 
£ 

Zr = 1 + Bi^ + + 0{u') , 

e £ 

£ £ 

Z^=l + CA + + 0{u^) , (43) 

Z^ = l + D^-L + D2-^ + /^3if^ + 0{u^) . 
2w£ 2w£ 2w£ 

Here, O(m^) is short-hand for the two-loop corrections which are of fourth order in 

the couplings u\\ and u±. Defining (o = i^di^ln Zo\hare = ~^{''^\\du^^ + u±duj^)zi^\ 

(y = — |(m||(9„|| -|- u_Lduj_)Y^-^^ as before, one obtains to one loop order 

C = -2(1 + w)ul , Cr = 3mJ - (5 + 2w)ul , 
(y = 3mJ - 3(7 + lOw + 4w^)ul , 

= (3 + w)ul - (5 + 9w + lOw^ + Aw^)ul , (44) 
(7 + 5w) ^ (49 + 65w + 26w;2) „ (10 + 2«;) 
^ lOw II lOw ^ lOw II 

We now consider the Gell-Mann-Low functions for the three dimensionless effective 
couplings u\\, u±, and w: 

2 + 8C-Cr+gCM-^(YT«;). 
/^.-'^^.^A.eH-| + ^C+^C.-^)n. (45) 



(3^ = K(9«w|bare = (^C - Cx + ^C)^ 



We seek a set of fixed points {w* ,u*^,u*j_) such that all three /^-functions vanish 
simultaneously. The eigenvalues of the linearized system, evaluated in the vicinity of 
each fixed point, provide information about its stability. If all of them are positive, the 
fixed point is infrared stable; otherwise, it is unstable in one or more directions. 

There is only one nontrivial, physical (i.e., real) infrared stable fixed point given by 



w* = 2^^-l + Oi£) , = ±y Mi±^ + 0(^3/2) ^ ^ Q 

For = 0, we find another root vu* = —1 — |"\/T5 < —1, but this solution must be 
excluded since it leads to a linearly unstable theory (see the discussion after Eq. ©). 
There are numerous fixed points with u*j_ ^ and u*^ 0, but these are all infrared 
unstable or nonphysical, in the sense that either m| or u*j_ are imaginary or w* < — 1. 
In contrast to earlier results by ^7], we do not find a physically viable fixed point with 
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■u^ 7^ 0. We have also carefully considered the possibility of "degenerate fixed points" 
with w* = or w* = oo; none of these are stable. 

Returning briefly to Eq. ()42|) . we note that the critical parameter t\\ appears only 
in the combination f|| = /x^^^^r|[, which is actually the appropriate scale-invariant form. 
This suggests a slight modification of the general renormalization scheme, Eq. p7|l . so 
that the renormalizations of the critical control parameters are now written as 

f II = Z-'/'Z;'/\Z^f\\ + Yt^) , = Z-W^ , (47) 

Their derivatives can be expressed in terms of the ^-functions of Eq. ()44|1 : 

-c , c 



'^'^.^lllbare = ( 2 ^ f " ~ ^^^^ ' 



.2 2 

«:9.r^lbare =Cr±- (48) 

One can see quite easily that = is an invariant subspace of the RG, since the 
coupling alone cannot generate any corrections to u*j_. Assuming that the stable fixed 
point, at higher orders of perturbation theory, remains characterized by w* 7^ 0, u| 7^ 0, 
while u*j_ = 0, we can establish the (exact) relations 

C = 0, lQ-2C = e, C:=Iq (49) 

which will lead to exponent identities, as we shall presently see. At the fixed point, the 
RGE takes the simplified form 

^4 - + ((iC - Of, - Qr.) I-] G.,^ = . (50) 

The two critical control parameters (f||,r_L) are diagonalized by introducing the new 
independent variables (o"||,r^), defined by 

^11 = ^11 + 2C*-C*/ ^ ■ ^^^^ 

This new critical control parameter cry takes over the role of ry in the mean-field 
approximation. The solution of the RGE ()50|) is given by 

= Fj, j^{{r, t}; r^; r^^/iy; 7, Ik) , (52) 

where I is an arbitrary scale, and -F/vtv results from Gj^fj after the variable 
transformation to {cr\\,Tj_). Combining this equation with the simple parallel and 
transverse scale invariances, we arrive at the scaling form 

G^,^({r, t}; r||, rx) = /V^F^^^({/i+Aa;||, /r^, IH}; r^/^^^a^^^l-^/^^T^) , (53) 

where fj is again given by Eq. while the other exponents take the values 

A r7 = r = 0, ^ = 4-r/ = 4, 

1M = 2-| + C, 1M = 2. (54) 
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The result to all orders, Eq. (|49|) . leads to an exact relation between A and u^: 

= r ■ (55) 

Once again, only a single exponent must be computed explicitly within the e-expansion, 
e.g., A. Our one-loop calculation yields 



. 23 + 6v 15 „ , 2\ /rn\ 

A = e + 0{e^), (56) 

and for the mixing of the critical control parameters 

^11 = ^ + f .../^ /T^. + 0{e))T^ . (57) 



Ti\\ V3(10 + V15) 



z/|l < - , A > . (5^ 



Let us briefly return to the assumption that the stable fixed point is characterized 
by = 0. Thanks to Eq. ()45|) and the associated linear stability matrix, we can 
determine the stability criterion, i.e., the correction-to-scaling exponent of m^, to all 
orders in perturbation theory. The condition for = to be stable is As < 3(* which 
generates bounds for the critical exponents A and , namely, 

1 , 6-d 
A > 

2 ' 3 

If the stability condition Ae < 3^* or equivalently, the second bound of Eq. (j3H|) . 
is violated at some dimension d, the fixed point structure of the theory changes 
fundamentally. Now 7^ is stable, and Eq. gives us the (exact) relations 

4e-3C; = 9C>0, Cx = ^(C + C) (59) 

instead of Eqs. (jl^jl . We lose one condition, and this leads to two independent critical 
exponents, as opposed to a single one. Following the same analysis as before, we arrive 
at the relations 

A = ^-^-r/, r]>0, l/iy^ = 2-r]. (60) 

The scaling behavior of the Green functions still follows Eq. but can no longer 
be related to A. 



4. The roughness exponents 

The roughness exponents of the surface are easily identified, once the scaling properties 
of the underlying field theory are known. For isotropic theories, they can be read off 
directly from the height-height correlation function, 

G2,o(r - r', t - t') = C{r -r',t- t') = {h{r, t)h{r' , t')) 

if its asymptotic scaling behavior of can be written in the form 

C(r,t) = |r|'^c(t/|rn (61) 



Novel surface universality classes with strong anisotropy 



19 



Here, x denotes the roughness exponent and z the dynamic exponent of the surface 
while c is a universal scaling function. Similarly, in Fourier space, the behavior of 
C{r,t) translates into 

C(q,t) = |qr('^+^^)c(t|qr) (62) 

with the same roughness exponent x- our case, the situation is slightly more subtle 
[SB], due to the presence of strong anisotropy First of all, surface fluctuations along 
the parallel and the transverse directions in real space need not be controlled by the 
same indices, leading us to define two exponents, x± and via 

C(x||,r^ = 0;t)= c,, (63) 

A similar situation should be expected in Fourier space, prompting us to define two 
additional exponents, x± and x\\, via 

5(0,qx,t) = |q±r^'^'*"^c^(t|q±r) 

C{g,,0,t) ^g,/^^^*"^c„(tgj/(^^^)) (64) 

Of course, all of these expressions are only meaningful if the four scaling functions 
c±, C||, Cj_ and c\\ approach finite and non-zero constants when their arguments vanish. 
Generically, the four roughness exponents take different numerical values; however, they 
can all be expressed in terms of the exponents A and rj, according to the identities 

X^ = ^[A-{d + A)-v] XII = ^[4-(rf + A)-r^]/(l + A) 

X± = l[^-d-r^] XII = ^[4 - + A) - ^ (65) 

The key observation is that xi. = X\\ = X± = X\\ only if the anisotropy exponent A 
vanishes. 

In the following, we explicitly compute the roughness exponents for the two cases 
(i) and (iii). For case (ii), the renormalization group gives us no such information, 
since that would require the existence of an infrared stable fixed point. If, in fact, 
the presence of a first order transition line were to be confirmed eventually, the whole 
concept of scaling exponents would be inapplicable here. 

Case (i) For this theory, characterized by positive t\\ while t± vanishes, we found that all 
exponents could be expressed in terms of rj = —^e + 0{e'^). Writing all four roughness 
exponents in terms of this single index, we arrive at expressions which are exact to all 
orders in e = 3 — d: 

l-d/3 



X± = l-d/3, x\\ 



3-d/3-r]' 



Xx = (4-d-,)/2. X1- ,(3,';;'_^) -^. (66) 

The mean-field values are easily recovered by setting e = 0. The physically most 
interesting case corresponds to a surface grown on a two-dimensional substrate, i.e.. 
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d = 2 and e = 1. For this situation, one obtains x± = 1/3 and x\\ = 1/9 + 0{e'^) while 
X|| = —2/9 + 0{e^) and x± = 4/3 + 0(£:^). Remarkably, the exponent x± = 1/3 is 
actually exact, at least to all orders in perturbation theory. 

Case (iii): At the multicritical point one has r] = 0, z = 4, and a nontrivial A. Hence, 
one obtains 



Xx = |4-d|/2 *,= (^-^ (67) 

All of them are negative near the upper critical dimension dc = Q. In order to access 
the physical {d = 2) situation, one has to set e = 4 here which gives a huge anisotropy 
exponent, A ~ 16.814, if one naively uses the one-loop results. While roughness 
exponents can in principle be calculated, it does not appear meaningful to compare 
them to experimental data. Yet, we emphasize again that the exponents calculated here 
do not agree with [T7j . 



5. Conclusions 



To summarize, we have analyzed the effect of strong anisotropies on the universal 
behavior of a surface grown under ideal MBE-type conditions. The anisotropy is 
generated by an incident particle beam which is tilted away from the normal. This 
arrangement manifestly breaks rotational invariance in the substrate plane. As a 
consequence, the effective surface tension becomes anisotropic, being characterized by 
two parameters, ry and t±, as opposed to a single one in the isotropic case. If both 
are positive, the surface is described by the Edwards- Wilkinson model. However, as 
temperature or other growth parameters are varied, the two control parameters need 
not vanish simultaneously, leading to different dynamic and spatial instabilities and 
ordered structures. Clearly, there are three distinct scenarios: (i) ry > while r_|_ goes 
to zero, (ii) ry —>■ while t± remains positive, and (iii) both parameters vanish. Each 
of these three cases is described by its own characteristic field theory with distinct 
upper critical dimensions: dc = 3, 7/2, and 6 for cases (i), (ii), and (iii), respectively. 
To display scale-invariant behavior in the long-time, large- distance limit, a field theory 
must possess one, or more, infrared stable fixed points. Our one-loop analysis reveals 
that only cases (i) and (iii) have this property. In contrast, case (ii) may correspond 
to a line of first-order phase transition whose characteristics lie outside the scope of 
our RG techniques. For the two scale-invariant theories, we find two distinct, novel 
surface universality classes, and derive the associated scaling behavior of the height- 
height correlation function. An intriguing feature of both universality classes is the 
emergence of a strong anisotropy exponent, which reflects different scaling behavior 
in the parallel and transverse directions. As an immediate consequence, we need to 
introduce four different roughness exponents, to characterize surface fluctuations in real 
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and momentum space, observed along the parallel vs the transverse directions. While 
all four are related by exponent identities, care must nevertheless be exercised when 
analyzing experimental data. 

We conclude with a few comments on the phase diagram. The first case clearly 
corresponds to a second order line, parameterized by ry. Assuming that the first 
order nature of the second case, parameterized by r^, can be confirmed, the third case 
T"!! = Tj_ = would turn out to be a critical endpoint. Since it requires the tuning of 
two parameters, it may be quite difficult to access in a typical experiment. In contrast, 
either the second or the first order lines should be more easily observable. From an RG 
perspective, even if the TySy /i term is absent initially, it is immediately generated under 
RG transformations, leading to a non-zero critical value of t\\; this is not the case for 
the TxV^h contribution. For this reason, we believe that the most physically relevant 
theory (beyond Edwards- Wilkinson behavior) is the one with ry > and t± = 0. Our 
key results for this model are the roughness exponents for real-space surface scans. For 
a two-dimensional surface embedded in a three-dimensional space, we find x± = 1/3 
for scans along the transverse direction, and x\\ = 1/9 + 0(5^) for parallel scans. There 
appears to be some experimental evidence for these exponents jHEI but more data are 
needed before this issue can be settled. 

We thank U.C. Tauber, R.K.P. Zia, J. Krug, A. Hartmann, and E. Yewande for 
helpful discussions. This work is partially supported by the NSF through DMR-0308548 
and DMR-0414122. GP acknowledges the Alexander von Humboldt foundation for their 
support. 

Appendix 

In the following, some technical details of the above calculations are presented. We prefer 
the time-momentum representation of the propagators and correlators, since it displays 
the causal structure of the Feynman diagrams most directly. In this representation, 
each Fi^i- (or self energy) diagram has a single h-\eg at the outgoing vertex with the 
largest time argument, and a single h-\eg at the incoming vertex with the smallest time 
argument. To find the UV divergent u; = contribution, we just integrate over the time 
difference of these two vertices. 

Appendix A. Case (i): t± ^ and r ^ 

To characterize the theory fully, it suffices to calculate the vertex function Fi i. The 
only one loop diagram which must be evaluated is 



using the diagrammatic building blocks associated with Eqs. and (fT^ . This graph 
represents the one loop contribution to the self-energy. After integrating over the time 




(A.l) 
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difference of the two vertices, it takes the form 

__ 2 f d'^p [-(qi)'p; + 2qi(q±-p±)pg + g|(pi ^^ 



\y J (2^)^ r(p) [r(p) + r(q - p)] . ^^'^^ 

Several finite terms have been dropped, following a dimensional regularization and 
minimal subtraction scheme. The function r(p) contains, up to a factor 7, the full 
momentum dependence of the inverse propagator 

r(p)=pi(pi + rx)+pfr|| . (A.3) 

Dimensional analysis immediately suggests that the terms in the bracket of 
Eq. ()A.2jl diverge logarithmically, except possibly the term 2q^(q^ ■ p±)p|, which is, 
at least naively, linearly divergent. The denominator in Eq. ()A.2|) is now expanded in 
small external momenta, and multiplied with the numerator, keeping only divergent 
contributions. After some straightforward algebra, we arrive at 



21) = -7A 



2 



i(qi)^/(2, 0, 2) + igJ/(0, 4, 2) + ;^(qi)^/(2, 4, 3) 



The integrals 7) are defined as 



(A.4) 



/(a,/5,7) 



with the proportionality indicating the degree of divergence based on simple power 

counting in d = dc = 3. Clearly, all integrals appearing in Eq. ()A.4|1 are only 

logarithmically divergent; the (naive) linear divergence vanishes by symmetry. Using 

the effective coupling u introduced in Eq. (j^U, the integrals are 

11^ 11^ iP' 

Ai/(2,0,2) = -, Ai/(0,4,2) = -r||, Ai/(2, 4, 3) = - , (A.6) 

and one finally arrives at Eq. ()2()|1 . 
Appendix B. Case (ii): ry and 7^ 

The diagrammatic representation of Fi 1, to one loop, is again given by Eq. dH}. 
However, the relevant interaction vertex here is Ay, so that the one loop integral reads 

7AM/(g;Pgfaii-Pii) V(p)|r(p) + r(q-p)| 

and the momentum dependence of the inverse bare propagator is controlled by 

r(p)=pj(/x||pj + r||) + r^pi. (B.2) 
As in the previous section, the integral ()B.1|) cannot be performed in closed form. Again, 



we expand the denominator of Eq. ()B.1|) in powers of qn and q^, noting that higher orders 
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become more and more convergent. Keeping only divergent terms, the singular part of 
Fi 1 takes the form 

(2^)d2r(p)2 (^-2) 



The leading singularity, i.e., the first integral above, is found to be quadratically 
divergent. It contributes to the renormalization of ry and can therefore be controlled by 
an additive renormalization. We shift the critical parameter by an amount f\\^c, defined 
through the vanishing of the parallel susceptibility, Eq. (jl5j) . Evaluating x\\ to first order 
yields 

xr-„-Ai/^4^ + 0(Aj). (B.4, 
The general renormalization scheme, Eq. (fTTj) . then allows us to determine 



he = ^ 

r||=0 



(27r)^2r2(p) 



0(A). (B.5) 



in one loop. Reparameterizing Fi^i in terms of the shifted ry subsequently removes the 
quadratic divergence and leaves us with logarithmic divergences only. Hence, the vertex 
function reads, to first order. 

El,! =i^ + 7(gy(/i||gy+r||) + rxql) + 7AjgJf|| J(6, 3) 

- ^7A|gJ [J{2, 2) - 7/i|| J(6, 3) + 4/iJj(10, 4) + ...] 

with the integrals J(q;,7) defined by 

= f -^^V oc A2<i-i+-47 (B.6) 
^ ' " J (27r)'^F(p)^ ^ ' 

Again, the proportionality indicates the degree of divergence based on d = d^ = 7 /2. 
Using the effective coupling u introduced in Eq. (jSHI), the integrals can be evaluated: 

Aj J(2, 0, 2) = , Aj J(6, 0, 3) = - , Af J(10, 0, 4) = /^^'^ , (B.7) 

which finally results in Eq. (jSTj)- 



Appendix C. Case (iii): r and t± ^ 

The calculation of the self energy follows the same standard methods as above. However, 
the technical details become more involved, for two reasons. First, the full momentum 
dependence of the inverse propagator, 

r(p) = pKpI + r^) + 2pipJ/ix + pjif^wpl + n) (c.i) 
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comes into play, with the two critical parameters, t± and ry , serving as infrared cutoffs. 
Second, both nonlinearities, parametrized by Ay and A^, are now marginal at the upper 
critical dimension dc = Q. 

As before, we start with the diagrammatic representation, Eq. ()A.1|) . Due to the 
presence of the full vertex, Eq. ()12|1 . the one-loop integral is given by 

f d'^p V{{q,p}) 

V (27r)'^r(p)[r(p) + r(q-p)] 

where V^({q, p}) summarizes the momenta contributed by the vertices: 

n{q,P})= -Ajg;pJ(g||-p||)' 

+ \\{ [2(q^ ■ px)p|| - q^pii - giipi] [q^pn - gup^] } 
+ A||Axg|| {q\\ - p\\) p\\ [2(q^ ■ p±)p\\ - 2qlp\\] 

As for case (ii), a shift of the critical parameter ry absorbs a quadratic divergence in 
Fi 1. Demanding that the inverse (parallel) susceptibility, Eq. (fT3j) vanishes at the true 
critical point, one finds in one loop order 

d'^p Pj 



n,c = All 



(27r)^2r(j9)^ 



(27r)^2r(p)^ 



= TX=0 



-Ai 

=rx=0 

Again, we retain only logarithmic divergences in Fi i if we reparameterize it in terms 
of a shifted ry ry — rn c. Expanding the denominator to O(q^) results in a sum of 
integrals of the form 

d'^p Pf(pi)^ /^ 

{2nY F(p)" 

Formally, these look like Eq. ()A.5|) . but it is essential to realize that F(p) takes the more 
complicated form ()C1|) . To evaluate the /(a,/3, 7), it is convenient to introduce a new 
integration variable, x = Q'||/|q±|, to replace gy. This has the significant advantage 
that all e-poles are already explicitly displayed, once the q^-integration has been 
performed. The remaining integrals over x can be evaluated at £ = and are perfectly 
straightforward. They take the form 

/oo 2a 
[l + 2wx^ + x^) 

with (n, a) e {(1, 0), (2, 0), (3, 0), (4, 0), (2, 1), (3, 1), (4, 1)}. Clearly, these K's lead to 
w-dependent coefficients for the different contributions to Fi 1. The result is summarized 
in Eq. ()42|) with the w-dependent functions 

A =2{l+w) , 

fii = -3, B2 = {5 + 2w), Bs = 3{7 + 2w{5 + 2w)) , 

= - (3 + w) , C2= {5 + 9w + lOw^ + Aw^) , (C.5) 

Di = -i(7 + 5w), 1^2 = 7(49 + 65^ + 26^2), D3 = I {10 + 2w) . 
5 5 5 

Tallying up all the contributions to Fii, and removing the divergences by 
renormalization, one arrives at the Z-factors listed in Eq. (j43|) . 
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